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Abstract. A constitutive model is developed to predict the viscoelastic response of polyimide 
resins that are used in high temperature applications. This model is based on a thermodynamic 
framework that uses the notion that the 'natural configuration' of a body evolves as the body 
undergoes a process and the evolution is determined by maximizing the rate of entropy production 
in general and the rate of dissipation within purely mechanical considerations. We constitutively 
prescribe forms for the specific Hclmholtz potential and the rate of dissipation (which is the 
product of density, temperature and the rate of entropy production), and the model is derived 
by maximizing the rate of dissipation with the constraint of incomprcssibility, and the reduced 
energy dissipation equation is also regarded as a constraint in that it is required to be met in 
every process that the body undergoes. The efficacy of the model is ascertained by comparing 
the predictions of the model with the experimental data for PMR-15 and HFPE-II-52 polyimide 
resins. 



1. Introduction 

Polyimides are well known to be extremely stable at high temperatures and also have a glass 
transition temperature that is greater than 300°C. Due to their good performance at high tem- 
perature ranges they are used by aircraft and automobile industries to fashion their products. 
They are also used in w afer fabrication due to th eir excellent high temperature resistance and 
adhesive properties (see Ghosh and Mittal ( 1996| )). The mechanical properties of polyimides 



and polyimide composites used in several applications especially in the aerospace industry are 
affected by high temperature, diffusion of moisture and subsequent oxidation. Hence, there 
is need for a good understanding of the various degradation mechanisms that are operational 
when such materials are subject to hostile environment. Recent experiment al evidenc e show s 



that the response of polyim i de re sins is solid-like viscoelastic response (see iBhargaval ( 120071 ). 



Falcone and Ruggles- WrennI (120091 )1 Moreover, the response of such bodies is non-linear (For a 



detailed description concerning how to differentiate between vis c oelastic solid-like and fluid-li ke 
response we refer the reader to Wineman and Rajagopal ( 2000l ). Malek and Rajagopal ( 2006 )). 



A thermodynamic framework which takes into account the viscoelastic solid-like response of the 
polyimides along with the various degradation processes needs to be developed. As a first step, 
our aim in this paper is to develop a model based on a thermodynamic framework that can predict 
the non-linear viscoelastic solid-like response of polyimides at various temperatures. This first 
step is non-trivial and presents interesting challenges. Subsequently, we shall extend this model 
to include degradation due to moisture diffusion, and chemical reactions, specifically oxidation. 
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While the thermal response of linear viscoelastic solids have been studied in great detail, 
there has been no systematic study of non-linear viscoelastic solids. Standard techniques like 
superposition that are valid in linear response are no longer vali d, thus making the s tudy much 
more complicated. A single integral model has been proposed by iPipkin and Roger 1 (ll968l ) who 
have assumed that a linear combination of responses to single step strain histories can be used 
as an approximation to the response to an arbitrary strain history. Unfortunately, such a model 
does not have a sound thermodynamic basis, and moreover the model is too general to be of use. 
Later on, iFung fll993h developed a quasi-linear viscoelastic model that has been shown to predict 
the behavior of several biological materials, though not adequately when the strains are large. 
This model by Fung can be shown to be a special case of the model by Pipkin and Rodgers. 
For further details on the various v iscoelastic mod e ls for solid s that hav e been reported in the 
literature, see the review articles by Drapaca et al.l ( 2007 ) and Wineman ( 2009 ). 

Most of the literature con cerning; the mod eli ng of the response o f poly imides use a vis- 
coelastic model proposed by Schapery ( 1969I ). Muliana and Sawant ( 2009 ) used Schapery's 
model and obtained material parameters for PMR-15 Q using the experiments carried out by 



Marais and Villoutreixl ( 1998 ). They used these material paramet ers in their 'micro - mecha nical' 



model to predict the behavior of Kevlar/PMR-15 composites. Ahci and Talrejal ( 2006 ) per- 
formed experiments on a composite made of graphite fiber in a HFPE-IlB They have extended 
the framework developed by Schapery to include a 'damage tensor' as an internal variable, and 
have also included anisotropy to model the composite behavior. Recently, Falcone and Ruggles- 
Wrenn performed experiments on PMR-15 at the service temperature of an aircraft, 288°C, and 
compared the pr edictions of Schapery's model for the problem of creep with experimental data. 
BhargavaJ ( 2007) has also used Schapery's model to predict the behavior of HFPE-II-52. More 



recently. iHall ( 20081 ) developed a thermodynamic framework for finite anisotropic viscoplastic 
models to study the response of polymers subject to extreme thermal environment. 

Given its extensive use some comments on Schapery's model are warranted: Schapery devel- 
oped a viscoe lastic model us ing linear phenomenological relations based on Onsager's reciprocity 
theorem (see Onsager ( 193ll )) which states t hat the fo r ces ar e linearly related to the fluxes near 
equilibrium (see for instance, equation 11 in Schaperyl ( 1964 )). Next, he introduced nonlinearity 
by assuming that the coefficient matrix relating the forces and fluxes depends on generalized 
coordinates and temperature. Furthermore, the free energy expression was obtained using a 
Taylor series expansion and by neglecting higher order terms. Thus, while the model might 
be able to describe slight deviation from linear response, one cannot expect it to be capable 
of describing truly non-linear response that is thermodynamically compatible. Thus, if one is 
interested in describing the non-linear response of viscoelastic solids that takes into account its 
thermodynamic effects, a different model is necessary. In this paper we develop a viscoelastic 
solid model based on a thermodynamic framework that can be used to describe the non-linear 
response exhibited by a class of polymers. The framework has been recently developed and 
is used to describe the response of bodies that produce entropy in a variety of ways. In or- 
der to derive meaningful physical models they require that amongst the class of processes that 
are possible the process which is actually taken by the body is one that maximizes the rate 
of entropy production. One can find details concerning this approach in the review article by 



PMR polymerization of monomer reactant. 
2 HFPE stands for hydrofluoropolyether. 
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Rajagopal and Srinivasa ( 2004al ). In this approach, one need not assume near-equilibrium be- 



havior and linear phenomenological relations between forces and fluxes, the approach is much 
more general. Also, one need not use a Taylor series e x pansio n of the free energy, and neglect 



higher order terms. Recently, iRajagopal and Srinivasa! (|2004bJ) have shown that if one uses an 



expression for entropy production which is quadratic in the fluxes, one can arrive at Onsager's 
relations upon maximizing the rate of entropy production along with appropriate constraints. 
As mentioned earlier such a thermodynamic framework has also been used to model various 
material responses such as viscoelastic solid-like and fluid-lik e behavior, traditional p l asticity , 



twinning, crystallization and so on (see the review article by IRajagopal and Srinivasal ( l2004al ) 
for the references and for the details of the framework). 

In this paper, a viscoelastic solid model is derived by assuming forms for the Helmholtz poten- 
tial and the rate of dissipation, and maximizing the rate of dissipation with incompressibility 
and the reduced energy dissipation equation as constraints. This model is shown to predict 
the viscoelastic response of polyimide resin. Experimental data for PMR-15 poly imide resin 



from Falcone and Ruggles- Wrenn ( 2009 ). and for HFPE-II-52 from Bhargava ( 20071 ) are used to 



evaluate the efficacy of the model. 

The current paper is organized as follows. In section (J2]), the kinematics that is required in 
this paper are documented. In sections (13.11) . ( 13. 2p . a viscoelastic solid model is developed using 
a thermodynamic framework. We show that the viscoelastic solid model that is developed is a 
generalization of the one- dimensional standard linear solid model in (13.31) . In section f!3.4[) the 
problem of uniaxial extension is set up using our model, and the creep solution obtained by using 
the model is compared with experimental data for PMR-15 and HFPE-II-52 polyimide resins in 
section (13.51) . We find that the theoretical predictions agree quite well with the experimental 
results. 

2. Preliminaries 

Let kr(B) and ft t (<B) denote the reference configuration and the current configuration, respec- 
tively. The motion \k r is defined as the one-one mapping that assigns to each point X G hr, a 
point x G hi t , at a time t, i.e., 

x = XKR (X,t). (2.1) 

The mapping Xk r (X, t) is assumed to be sufficiently smooth and invertible. Let k p ^) be the stress- 
free configuration instantaneously reached by the body upon removal of the external stimuli. We 
assume that the body can be instantaneously unloaded. We shall call this configuration as 
the natural configuration corresponding to K t . The natural configuration that underlies the 
current configuration depends on the process class that is admissible. Thus underlying natural 
configuration corresponding to isothermal and adiabatic processes could be different. Let F be 
gradient of motion x KiJ (X,t) (usually known as the deformation gradient), defined by 



9 Xk r 



(2.2) 



dx ' 

and let the left and right Cauchy-Green tensors be defined through 

B = FF T , C = F T F. (2.3) 



3 In case of isothermal processes, the rate of dissipation is the rate of conversion of mechanical working into 
heat (energy in thermal form), and in general it is the product of density, temperature and the rate of entropy 
production. 



Let F Kp{t) be the gradient of the mapping from k p (^ to Kt, and let G be defined by 

Similar to (12.3j) . we shall denote the left Cauchy-Green stretch tensors B G and B p (t) as 

B G :=GG T , B p{t) :=F Kpit) F T Kp(ty (2.5) 
We shall also define the velocity gradients 

L G :=GG~\ L = FF\ L p = F Kpit) F-^ (2.6) 
and their symmetric parts by 

Di = - [Li + Lj) , i = p(t), G or no subscript. (2.7) 

Also, we define the principal invariants through 

I fll = tr(B,), II Bi = ^{[tr(BO] 2 -tr(Bf)}, = det(B,) l = G,p(t), (2.8) 

where tr(.) is the trace operator for a second order tensor and det(.) is the determinant. Now, 
from (J23D 

F = F *pm G + F K P (t) G 

=> FF 1 = F Kp(t) GG-'F-l {t) + F Kp(t) G (2.9) 

=► L = L p[t) + F Kp{t) L G F- l m , 
where (.) is the material time derivative of the second order tensor. In addition, 



B m = F F T + FF T 



(2.10) 



- L p(t) B p(t) + B„(t)Xp( t ), 
and similarly 

Bg = LqBq + BqL g . (2.11) 
Hence, from (12.91) and (I2.10p . we have 

B p (i) = L% + B p{t) Lp {t) - F Kp{t) (L G + ^g) ^ t) i ( 2 - 12 ) 

and so 

B p(t) =-2F S(t) B G F^ w , (2.13) 

v v . 

where (.) is the usual Oldroyd derivative defined through A: = A — LA — AL . When one 

considers non-isothermal processes the local form of the second law of thermodynamics takes the 

following form: 

T ■ D — q%1) — qs6 — qh ' gT Q am = q9( := £ > 0, (2.14) 

a 

where T is the Cauchy stress, ip is the specific Helmholtz potential, q is the density, 9 is the 
temperature, s is the specific entropy, q h is the heat flux, ( is the rate of entropy production and 
£ is the rate of dissipation. 
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3. Constitutive assumptions and maximization of the rate of dissipation 



3.1. General results. We shall assume that the viscoelastic solid is isotropic and incompressible 
with the specific Helmholtz potential of the form 



i> = ip(B p(th B G ,9) = $(I B ,II B ,I Bg ,II Bq , 



(3.1) 



Since the elastic response is isotropic, without loss of generality, we choose K p (t) such that 



F = V 



0' 



(3.2) 



where V Kp{t) is the right stretch tensor in the polar decomposition of F Kp(t) . We shall also assume 
that the total rate of dissipation can be split additively as follows 



T ■ D - gif) - qsO = Cm > 0, 



q h ■ grad(fl) 

e 



£ c >o, 



(3.3) 



where £ m , £ c are the rates of mechanical dissipation (conversion of working into thermal energy) 
and dissipation due to heat conduction, respectively. Now, we constitutively choose 



q h = -k(8)grad(9), k{6)>0, 



(3.4) 



where k is the thermal conductivity, so that I3.3[ m is automatically satisfied. 
Next, 



dip 



^ | T 

dU m ^ B ^dII B , 



dip 



dip 



dll 



B 



B p (t) 



Pit) 



B 



p(t) 



01 



+ 1 



B G 



dip 



dip 



B G 



d\\ 



B 



G 



B G 



(3.5) 



and using fl2TT0l) . (jZHjl along with (E2D in (ESJ), we obtain 



ip = 2 
+ 2 



dip 



01 



B p(t) 



dip 



] 0ll 



B 



B p(t) 



pit) 



du~ 



01 



dip dip 
+ Ib g - 



Bg 



dll 



Bg 



dll 



Bg 



G 



(D-D 

dtp 



G 



■ D ° + M e - 



(3.6) 



Next, we shall assume the rate of mechanical dissipation to be of the form 



Cm — £m(#> 5p(t), D G ). 
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(3.7) 



On substituting (13.61) into (13.3^ )). we arrive at 

dip dip 



T-2q 

+ 2q 
-2q 



dl 



h 



B 



Pit) 



m dII B M 

dip 

8l B .„ ' ^ p(i) dll B ~ 



B P(t) + 2 P QII B l(t) 



D 



dtp 



I 



dip 



B n if\ — ~t~~~ — B 2 



B p(t) 



p(t) 



D 



G 



dip dip 



I 



di Ba ' ~ BG ~dii BG 
e 



- Q 



dip 2 



D 



G 



dip 



We shall set 



and define 



dip 

'He' 



T p (t) := 2g 



dip dip 



B 



p(t) 



dip 



B 



T G := 2q 



dip 



dlln p{t) 



+ 1 



d\ Bc ' ~ Ba ~dii BG 

Using f l3T9|) -f l3TTT]) in (E£g) . we obtain 

{T-T p(t) )-D+{T pit) -T G )-D G 
= €m(0, Bp(t), D G ). 
From constraint of incompressibility, we have 

tr{p) = tr(D p (t)) = tx(D G ) = 0. 
Since, RHS of (I3.12p does not depend on D, using A3 . 13j) . 

T = pl + T m , 

where p is the Lagrange multiplier due to the constraint of incompressibility, with 

(Tp(t) - To) ■ D G = £ m (0, Bp(t),D G ), 

which can be re-written as 

(T — T G ) ■ D G = £ m (9, B p (t), D G ), 

using fl3~T3|) and fl3~TD . 

Now, we shall maximize the rate of dissipation £ m by varying D G for fixed B p u\ . That is, we 
maximize the function^ 



dip 



B 



G 



(3.8) 



(3.9) 

(3.10) 
(3.11) 

(3.12) 

(3.13) 
(3.14) 
(3.15) 
(3.16) 



$ := U + Ai [£ m - (T - T G ) ■ D G ] + A 2 (J ■ D 



G 



(3.17) 



4 Though we only document that the first derivative is zero here, it can be shown that the extremum is a maximum. 
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where Ai, A2 are the Lagrange multipliers. By setting, d^/dD G = 0, we get 



T = T G + ^I 



Ax + 1\ d£ 



Ai J dD G 

We need to determine the Lagrange multipliers. On substituting (I3.18P into (I3.16p . we get 



(3.18) 



Ai + 1 
Ai 



dD c 



D 



G 



and so (13 . 1 8[) with (13. lip becomes 



T = 2q 



dip dip 
+ ls G - 



01 



B G 



dll 



B 



G 



B G 



dll Bo 



f£?- ■ D r I dD 



A/, 



(3.19) 



(3.20) 



where A := ^Ms the Lagrange multiplier due to the constraint of incompressibility. 
Finally, the constitutive relations for the viscoelastic solid are given by 



T = pi + 2g 
T = XI + 2g 



dip 



d A 1 T 

dl B m Bp(t) dll B 

dip dip 
+ Ib g - 



dip 2 



P (i) 



01 



Bo 



-k(9)grad{9), s 



dll 

dip 



B 



dip 



G 



B G 



dll 



B 



Be 



G 



din 

gas- • Dr ) 9D 



3.2. Specific case. Specifically, we choose the specific Helmholtz potential as 



iP = A s +(B s +c s 2 



9 s )-4(9-9 s ) 2 -cf > 9]n[- f )^ 



9 \ . A*G0 f^GlO 



2q9 s 



(Ib g -3)4 



(3.21a) 
(3.21b) 
(3.21c) 

(Ie p(t) -3), 
(3.22) 

where fi G , \i v are elastic constants, 9 S is a reference temperature for the viscoelastic solid, and 
the rate of dissipation as 

U = V (9) (D G ■ B p{t) D G ) , (3.23) 

where rj is the viscosity. 
Now, 

dip 

S = ~d9 

= -(B s + 4) + c{ (9 - 9 S ) + c s 2 \n (p) + 4 



2q9 8 



SI (Ibg " 3)+ 5: (Ib ^" 3) - 



(3.24) 



The internal energy e is given by 
e = Pp + 9s 



A' - B°9 S +4(9- 9 S ) + 91) + gL(L3 G - 3) + ^(I Bp(t) - 3). 



2q8 8 



and the specific heat capacity C v is 



Oc 



^ v ~ ~d~9 ~ C ^ C&2 ' 



(3.25) 



(3.26) 



Also, (E£2Iaj) . ( I3.21bp reduce to 

T = pi + p,pB p (t), 



XI + fi G B G + | {B p{t) D G + D G B p{t) ) 



where /i p 



fJ.pQ-fl.pl 8 - 



^. From OMB 



(p - A)/ + p.pBp(t) = [i G B G + - (B p{t) D G + D G B p{t) ) 



(3.27a) 
(3.27b) 

(3.28) 



and so by pre-multiplying the above equation by and taking the trace, we get 

pGtT(B~^B G ) - 3/Zp 



(p-A) 



(3.29) 



Using ( I3.29P in (13. 28 p . we arrive at the following equation that holds: 



fM G tT(B^ t) B 



*G) 



3fi } 



which can be re-written as 
'p G tr(B^ t) B G )-3fl 



I + P>pB p (t) = VgBq + - [BptyDa + D G B p ^) 



(3.30) 



VpBp(t) = fi G B G - - I V p (t) B p{t) V^ 1 + V^ 1 B p (t) V 



p(t) 



(3.31) 



where we have used ( 12.131) and ( 13.301) . Thus, with the current choice of the specific Helmholtz 
potential and the rate of dissipation, we arrive at the following constitutive equations: 



T = pl + fJLpB p (t), 



(3.32) 



where the evolution of the natural configuration is given by (13 .31 p. Also, note that the above 
model reduces to the generalized Maxwell fluid model derived by lRajagopal and Srinivasa (l2000h 
when p, G = 0. This is interesting, but not totally surprising, that we obtain a fluid model by 
eliminating a energy storage mechanism. In the corresponding one dimensional model this is 
tantamount to a spring being removed. 



3.3. Relationship to the Standard Linear Solid. Now, ( I3.32p . ( I3.3ip can be re- written as 

T = (p + fi P )I + fi P (B pit) - /), (3.33a) 



VGte{B p ( t) B G ) - 3/i f 



p(t)J 



I + pp(Bp {t) - I) — p, G {B G - I) 



V p( t) B p{t) 



>(t) 



Pit) 



v 
B 



P (t) v p(t) 



(3.33b) 



If Aj (i = G,p) is the one- dimensional stretch and = InAj (i 
when one is restricted to one- dimension, (I3.33P reduces to 



a 



fi p (X 



i: 



2 ^G 



G,p) is the logarithmic strain, 

(3.34a) 
(3.34b) 



where a is the one dimensional stress. Equation (13.341) under the assumption that <C 1 (i 
G,p) reduces to 



2fx p e p = 2ji G e G + rje G , 



(3.35a) 
(3.35b) 



where fj = rjX 2 is the stretch dependent viscosity. Equation (I3.35P can also be obtained by using 
a Kelvin- Voigt element (with spring constant 2jl G , viscosity of r/) and a spring (of spring con- 
stant 2p,p) in series, which is the spring-dashpot analogy for the standard linear solid. However, 
in this model the viscosity is stretch dependent and hence the model is a generalization of the 
classical standard linear solid as the viscosity in the standard linear solid model is assumed to 
be a constant. Hence, the viscoelastic solid model given by (13.321) . (I3.3ip is a three-dimensional 
generalization of the standard linear solid. Of cours e, there can be infinity of t hree dimen- 



Karra and Rajagopal ( 20091 )). Recently, 



sional generalizations o f a one dimensional model [see 

Kannan and Rajagopall (120041 ) have also derived a three-dimensional viscoelastic solid model, 
that is different from the model developed in this paper, that also reduces to the standard linear 
solid. That more than one, in fact, infinity of generalizations are possible is akin to the situation 
in elementary mathematics and stems from the fact that infinity of three dimensional functions 
can have the same one dimensional projection. In fact, even when one considers the thermody- 
namical formulation that is used in the paper, using different forms for the specific Helmholtz 
potential and the rate of dissipation, and by maximizing the rate of dissipation with the nec- 
essary cons traints more than one three- dimensional model reduces to the same one- dimensional 
model 



see 



Karra and Rajagopall (120091 ) for details of an example 



3.4. Application of the model. Let us study the uniaxial extension, given by 



x = X(t)X, y 



1 



=y, z 



i 



within the context of this model. The velocity gradient is given by 

L = diag < ——r, r 

& [A' 2A' 2A 

We shall assume that the stretch B p m is given by 

B K0 =dtag{fl,J=,-i r 



So, 



B p (t) = diag < B 



B 



' 2B 3 / 2 ' 2B 3 / 2 J 



(3.36) 



(3.37) 



(3.38) 



(3.39) 



V 

B 



p(ty 



diag < B 



2B\ 



B 



X 



B 



X 



V, 



A ' 253/2 Av /£' 2B 3/2 

= diag < \f~B. 



51/4' 51/4 



and 



G 



-diag 



5 



2A A B A 
T' A ~ 2B' A 



2£? 



Also, 



G 



K p(t) 



diag 



A B 1 ' 4 B 1 ' 4 



y/5' VA ' Vx 



which yields 



and 



Br 



A 2 



diag < — 



Substituting 1EE5B)) . ( ggpj) , (pj45j) into ( l3~30l) 



£' A ' A 
X 2 B B 



5 
2" 



BA /^ G A 2 // p 

X 4 



77 .B 77 

which can be re-written in the following form: 



A = A< 



I 2B 



1 
~rjB 



A 2 



^(A 3 + 2B 3 ) -3^AB 2 
XB (1 + 2B 3 / 2 ) 



/2 G (A 3 + 2B 3 ) - 3/i pJ B 2 A 



BX(1 + 2B 3 / 2 ) 



(3.40) 
(3.41) 

(3.42) 



(3.43) 



(3.44) 



(3.45) 



(3.46) 



(3.47) 



Now, from (13.381) . (13.321) . and using the fact that lateral surfaces are traction free, we conclude 
that 



(3.48) 



We shall also use logarithmic strain (or true strain) e = InA as our strain measure in what follows. 

3.5. Comparison with experimental creep data. For the loading process, with known con- 
stant applied stress Tu and material properties, ( 13.481) was first solved for B(t). Then, (I3.46H was 
solved with the initial condition A(0) = a/ B(0). For the unloading process, Tu was set to zero 
and B(t) was evaluated using (13.481) . Then, using A(i+) = as the initial condition (where t u 

is the time when unloading starts), X(t) during the unloading process is evaluated using (I3.46H . 
All the ODEs were solved in MATLAB using the ode45 solver. 
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In order to obtain the material parameters for a given set of experimental creep data, f minsearch 
function in MATLAB (which uses Nedler-Mead simplex method) was used to minimize the error 
defined by 



7 j y^theoload ^-expload) , /-, \ / 7 j (^-theo, unload ^-exp, unload) / Ar ,\ 

error = wx\l — -| + (1 - w) x W — ^ > ( 3 - 49 ) 

/ j \£-exp,load) y / y y^exp, unload) 

where Stheo denotes the theoretical strain values, e exp denotes the experimental strain values, 
the suffixes load, unload denote the values during loading and unloading processes respectively, 
w is a weight. The material parameters for the model was o btained for HFPE-II-52 polyimide 



resin using the experimental creep data from iBhargaval (120071 ) at different temperatures (285°C 



300°C, 315°C and 330°C). To determine the efficacy of the model the following process was 
followed. At 285°C, the experimental data values for the loading of 0.45 UTS were used to 
obtain the material parameters by minimizing the error (I3.49p . Then, these material parameters 
were used for the model prediction at the other loadings of 0.30 UTS and 0.15 UTS. The loading 
values corresponding to the sets of experimental data which were used to obtain the material 
parameters at the other temperatures are shown in table ([I]). Similar to the process described 
above for 285°C, the material parameters shown in table ([T]) were used to predict the creep at 
other temperatures. The model predictions compare well with the experimental data as shown 
in figures ©, ©. 

Next, the creep solution that stems fro m our model is compared to the experimental creep 



data of iFalcone and Ruggles-Wrennl ( 120091 1 for PMR-15 resin at 288°C in figure flU). The best- 
fit values of the parameters for were found to be fie = 4.42 x 10 8 Pa, ]x v = 3.76 x 10 s Pa, 
n = 6.22 x 10 12 Pa.s. A weight of w — 0.75 was used since there are fewer data points for 
the unloading process. As it can be seen from figure fll]), our model shows a good fit with the 
experiment. However, there is no additional experiment with which the predictive capability of 
the model can be tested. 

It is seen in the experiments that at a loading close to the failure values, the experimental 
data shows permanent set in the body, and there seems to be 'yielding'. Our model being a 
viscoelastic solid model it cannot predict such a permanent set. Thus, the model should be 
generalized to take into account the inelastic response of the polymer, but this is a daunting 
problem that requires a careful and separate study. 

In conclusion, a viscoelastic model has been developed which predicts the behavior of polyimide 
resins, that takes into account the thermal response, quite well. Our work can be extended to 
include anisotropy, and inelasticity to predict the response of polyimide composites and one can 
include various degradation mechanisms as well. 
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Table 1. Table showing values for the ultimate tensile strength (UTS) and various 
material parameters pc, 77). The table also shows the loading value data set 
that was used to obtain the optimum set of material parameters. 



reference configuration current configuration 




natural configuration 
Figure 1. Illustration of various configurations of the body. 
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Figure 2. Com parison of the model predictions with experimental creep data of 



Bhargaval (120071 ) at different loadings. The polyimide in this case is HFPE-II-52 at 
285°C and 300°C The parameters chosen and the values for the ultimate tensile 
strength (UTS) are shown in table ([I]). A weight of w = 0.5 was used for these 
two cases to obtain the optimum set of parameters. 



14 



2.5 



0.30 UTS (Model) 

0.15 UTS (Model) 

* 0.30 UTS (Experiment) 
a 0.15 UTS (Experiment) 



a3 



, *#*****************' 



0.5 



'aaaa a a a A A AAAa a AA A Aa a a a aa^ * 



V^Aaaaaaaaaa a,a/Yi^ m 



0.5 1 1.5 2 

time (s) 

(A) 315°C 



2.5 



x 10 



2.5 



1.5 



0.5 







0.20 UTS (Model) 

0.15 UTS (Model) 

* 0.20 UTS (Experiment) 
a 0.15 UTS (Experiment) 



********************** 



^aaaaaASASAaaaaaSaAaaaaaa 



0.5 1 1.5 2 

time (s) 

(b) 330°C 



2.5 



x 10 



Figure 3. Com parison of the model predictions with experimental creep data of 



Bhargaval (120071 ) at different loadings. The polyimide in this case is HFPE-II-52 at 
315°C and 330°C. The parameters chosen and the values for the ultimate tensile 
strength (UTS) are shown in table (fl|). A weight of w — 0.75 was used for these 
two cases to obtain the optimum set of parameters. 
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Figure 4. Comparison o f the model with experimental creep data of 
Falcone and Ruggles-Wrenn for a loading of 10 MPa. The polyimide in 

this case is PMR-15 at a temperature of 288°C. The parameter values used were 
p G = 4.42 x 10 8 Pa, fL p = 3.76 x 10 8 Pa, r) = 6.22 x 10 12 Pa. 
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